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INTRACLUSTER STARS AND THE CHEMICAL ENRICHMENT OF THE INTRACLUSTER MEDIUM 
Dennis Zaritsky\ Anthony H. Gonzalez^"', and Ann I. Zabludoff^ 

Accented for Publication in the Astrophysical Journal Letters 

ABSTRACT 

We explore the contribution of intracluster stars (ICS) to the chemical enrichment history of the 
intracluster medium (ICM) . In contrast to scenarios in which all the metals originate in cluster galaxies 
and are then transported into the ICM, intracluster stars enrich the ICM in situ, thereby contributing 
100% of their supernovae ejecta directly into the ICM. Modeling the ICS as an ancient, single burst 
stellar population with a normal initial mass function, we generate ICM iron abundances in the 
range of the observed values of several tenths solar. Large observational and theoretical uncertainties 
preclude us from concluding that the intracluster stars are the primary contributor of metals to the 
ICM in general. However, for the two clusters in our sample, and one from the literature, for which 
all of the required observational constraints exist, we are able to reproduce between half and all of 
the ICM iron with this simple model. Due to the ubiquity of intracluster stars in clusters and their 
direct connection with the ICM, we conclude that all models of the chemical enrichment history and 
energy budget of the ICM should account for the impact of the ICS. 
Subject headings: cosmology: observation — galaxies:clusters:general — supernovae 



1. INTRODUCTION 

The chemical enrichment history of the intraclus- 
ter medium (ICM) is arguably the least under- 
stood aspect of galaxy cluster evolution. Particu- 
larly puzzling are the relatively hi gh Fe abundance 
(~ 0.3 Fe^ : lEdsre fc S tewart 199 11) and moderatel 
large [a/Fe1 JRenzini et aL.1993;.Mushotzkv et alJll99 
iTamura et al.l l2004() . The former has been inter- 
preted as evidence for large gas outflows from galaxies 
l)R,enzini~ t al. 1993; Renzini 1997), while the latter, in 
combination with constraints from cluster galaxies, is 
interpreted as evidence for a non-standard initial mass 
function ( Mus hotzkv et all llQQfil: iPortinari et all 12004 
although see iPipino et all l)2002|) for a contrasting result) . 
Both of these inferences have a wide ranging impact on a 
variety of astronomical topics. Here we present evidence 
that the first argument should be significantly revised, 
and that the second argument should be revisited. 

Our proposal is straightforward. The intracluster stars 
(ICS) that we ha ve found without exception in a sam- 
ple of 24 clusters l)Gonzalez. Zabludoff. &: Zari tskv""2004l 
hereafter Paper I) must chemically pollute the ICM via 
supernov ae. Such intergalacti c SNe have already been 
observed l)Cal-Yam et al.ll2003D . We evaluate the degree 
to which the ICS plays a role with the simplest model- 
ing we can envision. As the vast literature on this topic 
demonstrates, there are many variables and unknowns 
in this type of model. We will adopt plausible, often 
typical, values of the relevant variables and demonstrate 
that the intracluster stars make a significant, and per- 
haps dominant, contribution to the enrichment history 
of the ICM. This study is not intended to be a complete 
or comprehensive treatment of the ICS enrichment of the 
ICM, but it illustrates the likely importance of the ICS 
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in the chemical abundance budget of galaxy clusters. 

2. OBSERVATIONAL CONSTRAINTS ON THE 
INTRACLUSTER COMPONENTS 

In Paper I we present observations of 24 clusters that 
span a range of velocity dispersions and Bautz-Morgan 
types (Bautz & Morgan 1970). The sample consists of 
nearby clusters (0.03 < z < 0.13) that contain a domi- 
nant brightest cluster galaxy (BCG) with a major axis 
position angle that lies within 45 degrees of the east- west 
axis (our drift scan direction). We present details of these 
unique data and our reduction procedure in Paper I. 

From modeling the two dimensional optical surface 
brightness distribution, we conclude in Paper I that the 
BCG light can be divided into two components. The 
inner component is similar to normal giant ellipticals in 
extent, luminosity, and scaling relations. The outer com- 
ponent is closely related in scale and ellipticity to the 
distribution of galaxies in the cluster. We associate the 
inner component with the BCG and the outer compo- 
nent with the ICS. In Paper I, we determine the lumi- 
nosity of the outer component, which we will use here as 
a measure of the intracluster stellar content. Standard 
measurements of BCG light underestimate the contribu- 
tion of this low surface brightness outer component (for 
example, assuming a single r^'^ surface brightness profile 
to evaluate the total magnitude leads to a 50% underesti- 
mate of the total luminosity of the BCG plus ICS; Paper 

I). 

For constraints on the properties of the ICM, we have 
searched the literature for X-ray gas masses ("Reiprichl 
HOOI) and chemical abundances (White 2000). There 
are only two clusters in our sample for which such data 
exist (Table QJ. We have adopted the gas mass inside 
of r^oo, the radius at which the mean mass density ex- 
ceeds 500 times the cosmological value, as the measure 
of the gas mass. We measure the ICS luminosity di- 
rectly only within 300 kpc, but our model fit provides an 
estimate of the total luminosity. The mismatched aper- 
^zal&iH(i<#@Bji^rtie«ip?<4i8 sion with which comparisons can be 
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made, although even if the apertures are matched, the 
comparison is comphcated by flows of material across the 
aperture during the cluster's evolution. The gas masses 
have been rescaled to our adopted cosmology {Hq = 70 
km/sec/Mpc, ftm — 0.3 and J^a — 0.7), and the ob- 
served abundances have been rescaled due to the differ- 
ence between the older valu e of the photospheric sola r 
abundance by mass (0.0026. lAnders fc Grevessgll989|) . 
and the more recent on e that we adopt here (0.0017, 
iGrevesse fc Sau var'1999"). This conversion increases the 
canonical Fe cluster abundance from 0.3 to 0.46 Fe©. 

3. MODELS AND RESULTS 

This type of study usually reli es on analytic mod- 
els of the enrichment hist ory dR enzini ct al." 1993; 
Brighenti fc Mathew^ 119981: iMaoz fc Gal- Yam 2004 : 
Portinari et alJ 120041) . but we choose to model the 



chemical evolution wi th the publicly available modeling 
package, PEGASE.2 IjFico fc Rocca-Volmcrangc 1997!), 
which provides estimates of the Type I and II super- 
novae rates as a function of age for a user-selected star 
formation history and initial mass function. We use these 
models because they also provide luminosities and colors 
of the resulting stellar population, which we need to scale 
our predictions for the observed ICS population. 

To mitigate ambiguity in our results, we adopt a very 
simple star formation history (a single burst of star for- 
mation that peaks at z > 10 and is finished by z = 4), 
an initial chemical abundance by mass fraction of 0.004 
for the ISM (although adopting values of 0.0001 and 0.08 
does not significantly affect the results), a standard IMF 
lIMiller fc Scalolll979l) . and all the default PEGASE pa- 
rameters (such as a binary fraction of 0.5). Changing 
the star formation history so that the peak star forma- 
tion rate is at z — 3.5 and the star formation does not 
cease entirely until z ^ 2 decreases the predicted abun- 
dances by ~ 10%. Varying the adopted initial chemical 
abundance only affects aspects of the stellar evolutionary 
modeling, such as the predicted stellar colors. We always 
assume that the ICM abundance is zero until we add the 
contribution from the ICS. 

By integrating the rate of SNe predicted from the PE- 
GASE models (Figure^ and adopting Fe yields for SNe, 
we obtain the total Fe mass ejected over time. For the 
yields we assume that a Type I ejects O.7M0 of Fe per 
supernova and that a Type II ejects O.O5M0 of Fe per 
supernova. This yield for a Type I is t he canonical value 
adopted by various other studies (see IMaoz fc Gal-YamI 
l2004HPortinari et alJ2004r. while that for Type IPs is the 
mean value fElmhamdi . Chugai. fc Danz iger"2003^ over 
the observed range, 0.0016 to 0.26 Mq (Hamuv 2003). 
Because Type I's dominate the enrichment in our mod- 
els, the results are fairly robust to decreasing the Type 
II yield from our adopted value. Increasing the Type 
II yield leads to more pronounced changes. However, in- 
creasing the Type II yield obviously produces a larger to- 
tal Fe mass and therefore is not a conservative approach. 

This discussion sidesteps the difhculties others have en- 
countered in re producing [g/Fe] or spec ific element ratios 
(jRenzini et al.lll993t iMushotzkv et al]ll996|) . However, 
any solution that has been proposed previously, such as 
non-standard initial mass functions, will also work with 
our ICS scenario, with the added flexibility that now one 
does not need to fit both the ICM and galaxy abundance 
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Fig. 1. — SNe rates for Type I (solid line) and II (dashed line) over 
time from the PEGASE models. The units are SNu (SNe per 10^° 
Ls per century). T he horizontal dotted line represent s the SNe 
rate observed locally ( Cap pellaro. Evans. &: 'LurattdlOS g) and the 
vertical dotted line represents the current time. For the adopted 
parameters, the PEGASE models produce a slightly conservative 
(lower) value of the SNe rate at the current time than is observed. 



patterns with the same model. Therefore, although we 
do not pursue the [a/Fe] question here, we conclude that 
it is no more difficult a problem than in the standard 
models. Future models should be able to take advantage 
of the wealth of information becoming a vailable on the 
spati al gradients of these element ratios IjTamura et alJ 
I2OOI . 

We present two different calculations that provide 
quantitative estimates of the importance of the ICS. 
First, we model clusters as a general population, adopt- 
ing typical ratios of such quantities as gas mass to stel- 
lar mass. Second, we model the only three clusters for 
which all of the required quantities are measured. In 
both approaches, we combine the PEGASE supernovae 
rates with the yields described above to calculate the 
Fe mass that 1 M0 of initial material produces after 
14 Gyr (this 1 Mq of material results in 0.59 M© of 
stellar mass, including remnants, at the current time). 
In the first type of calculation, we utilize the ratios of 
gas mass to galaxy mass in stars {Mgas/Mqal) and of 
ICS to total cluster luminosity (Ljcs/Ltotal), and as- 
sume that that the ICS and galactic stars have the same 
M/L ratio, to calculate the ratio of ICM gas mass to 
ICS mass. We then use this ratio, and the calculated Fe 
abundance per solar mass of ICS, to calculate an ICM Fe 
abundance. Values f or MgAR/MgAr, vary between ~ 2 
to 10 l|Arnaud et all 119921: iFukugita. Hogan. fc PcebQ 
Il998(l . although a combination of rec ent opt ical work 
with the Sloan Digital Sky Survey fBcU ct al. 200|) 
and a large ROSAT sample llEttori fc Fabian.1999 ) fa- 
vors the upper end of that range IjMushotzkvl 12004). 
In the limited observational work to date, values for 
Ljc rILtotaj, range from 0.1 (Peldmcicr ct al. 2004) to 
0.2 IjGonzalez et al.! '2000. adjusted for the division into 
two compon ents of Abell 1651 from Paper I). If, as mod- 
els predict IjSommer-Larsen. Romeo, fc Portinaril l2004t 



Chemical Enrichment of the Intracluster Medium 




Fig. 2. — Parameter choices that produce ICM metallicities of 0.1 
and 0.46 FeQ are shown by the Unes (an abund ance of 0.4 6 F ep, 
corre sponds to the canonical ^^ 0.3 abundance (Edge fc Stcw ara 
11991]) for our adopted photospheric Fe solar abundance). The ver- 
tical axis represents the mass ratio of the ICM to stars within 
galaxies. The horizontal axis represents the fraction of all cluster 
stars that are in the ICS. Although producing an ICM metallically 
of 46% solar strains the observational constraints on both axes (see 
text), producing a significant contribution, ~ 0.1 to 0.2 Fcq, is well 
within the observational limits. Given the modeling and observa- 
tional uncertainties, we cannot exclude the possibility that the bulk 
of the ICM Fe comes from the ICS. 



iMurante et al.ll200^ . the stars in the ICS are on average 
older than those in the galaxies, then our assumption of 
equal M/L for the ICS and galactic stars leads to an un- 
derestimate of the relative amount of mass in the ICS. 
This is therefore a conservative assumption in our deter- 
mination of the produced Fe. Instead of adopting specific 
values for the various ratios, we provide results over the 
range of plausible values (Figure [SJ- 

Figure 121 demonstrates that for the ICS to chemically 
enrich the ICM to cluster-like values requires both a 
fairly high Ljcs/Ltotal and a low Mqas/Mgal rel- 
ative to the ranges of values explored. However, to 
reach a level where the ICS contributes significantly 
to the ICM (~ 0.1 — 0.2 Fe©) requires values of the 
two scaling r atios that lie well within the observa- 
tional raneres llArnaud et al.l l l992t iGonzalez et al J l200(l 
iFeldmeier et al.ll2004j) . Furthermore, because M/L is 
likely to be larger for the ICS than for the galaxies (see 
above), our scaling may underestimate the contribution 
from the ICS by up to a factor of two. 

The second type of calculation provides a more direct 
comparison in the few clusters for which the relevant data 
are available. For two of our clusters, we find correspond- 
ing Mqas and metallicity measurements in the literature 
(Table ^. We scale the Fe mass, as calculated from the 
PEGASE models and our adopted yields, using the M/L 
from those models and the total luminosity of the outer 
component (Paper I) . The ICM without the contribution 
of the ICS is assumed to have zero metallicity. The cal- 
culated metallicity resulting from the ICS contribution 
is presented in Tabled for comparison. Mies represents 
the stellar mass, including remnants, in the intracluster 



volume. 

It is unfortunate that only two of our 24 clusters have 
published values of Mqas and Fe abundance. However, 
for those two clusters, we find that the ICS can account 
for a large fraction of the metals in the ICM, and we 
reproduce the relative rankings of metallicity for the two 
clusters. There are numerous adjustments that can be 
made due to uncertainties in the luminosity of the ICS, 
the region over which Mqas and the ICS are measured, 
the adopted yields, the modeling of the SNe rates, the 
star formation history, initial mass function, and so forth. 
In TableHl we also inclu de results for a cluster observed by 
iFbldmeier et al.l (|2004D . Despite some differences in their 
definition and measurement of the ICS, our model again 
reproduces a sizable fraction of the ICM iron abundance 
in this cluster. The purpose of this exercise is primarily 
to demonstrate that for plausible values of the various 
parameters, the ICS can contribute a significant amount 
of iron to the ICM. A similar conclusion has been reached 
using s caling relationships and infrared photometry by 
iLinfc Mohr (2004) . 

A complete model of cluster enrichment must include 
the contribution from all stars, whether currently in the 
ICS or in galaxies, and must account for the metals both 
in the ICM and in galaxies. Our approach has ignored 
1) the contribution of metals to the ICM from stars cur- 
rently in galaxies, because we are interested in deter- 
mining the contribution from the newly measured ICS 
component, and 2) the metals currently locked in galax- 
ies, again because we are interested in the connection 
between the ICS and ICM. In effect, we have treated 
the current cluster galaxies as closed box systems and 
set them aside. This approach produces a conservative 
estimate of the metal enrichment of the ICM. 

It is important to discuss one possible refine- 
ment of our chemical enrichment model. We 
make no direct assumption about the origin of the 
ICS, which, in simulations, arises from the tidal 
disruption of cluster galaxies early in the clus- 
ter's history fSommer-Larsen. Romeo, fc Porti narl2003: 
Murante et al. 2004; Willman et al. 2004). If these simu- 
lations are correct, then the metals in the ICM are some 
combination of the metals already present in those first 
cluster galaxies and the subsequent evolution of the stars 
stripped from those galaxies. So far in this paper, we 
have discussed only the latter issue. 

We now consider the former issue by asking what frac- 
tion of the metals in those ICS progenitor galaxies enters 
the ICM. If the parent galaxies of the ICS were com- 
pletely disrupted, then the answer is 100% and our model 
implicitly includes this possibility. If some of the parent 
galaxies survive to the current time, then our calculation 
overestimates the ICS metal contribution to the ICM by 
whatever amount of metals is locked in these galaxies. 
Without knowing when the ICS formed and its detailed 
star formation history, we cannot estimate the magnitude 
of this effect, but it may be significant. For example, for 
a star formation model that peaks sharply at ^ ~ 3, an 
age compa rable to the mean age of the ICS in the siin - 
ulation bv iSommer-Larsen. Romeo, fc Portinaril l)2004ll . 
about 25% of the Fe forms prior to the mean age of the 
ICS stars. If the stars are stripped from these galaxies 
shortly thereafter, the galaxies might retain as much as 
25% of all the metals that we attribute to the ICS. Nev- 
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TABLE 1 

Cluster Parameters and Estimates of the Intracluster Fe Abundance 



Cluster 



^BCG 



Mgas, 500 



rsoo 
(Mpc) 



Mr 



Ml 



Mies 
(IQi^Afo) 



Observed Fe 
(solar units) 



Predicted Fe 
(solar units) 



Abell 1651 
Abell 1914° 
Abell 3112 



0.0853 
0.1712 
0.0759 



6.5 
12.9 

4.5 



1.73 
2.46 
1.53 



-24.95 



-25.74 



-24.94 



0.9 
1.8 
1.9 



0.44 
0.37 
0.68 



0.22 
0.22 
0.66 



^This cluster is drawn from iFeldmeier et alJ 120041^ . Their measurement of the ICS is somewhat different than 
ours, partly because they observe clusters without a dominant brightest member. Wc adopt the ICS luminosity 
corresponding to their mask threshold, which best approximates our approach. Their measurement is in the V 
band, which further complicates the comparison. 



ertheless, even if our calculation overestimates the ICS 
metal contribution to the ICM by as much as 25%, the 
ICS still contributes between 35% and 75% of the metals 
in the three clusters that we model in detail (Table 1). 

4. CONCLUSIONS 

We calculate that the chemical enrichment due to in- 
tracluster stars (ICS), a ubiquitous population in galaxy 
clusters (Paper I), is a significant component of the 
chemical enrichment history of the intracluster medium 
(ICM). Furthermore, our model suggests that for at least 
some clusters the ICS could contribute the bulk of the 
metals in the ICM. Unlike models of ICM enrichment 
in which the metals originated in the cluster members 
we see today, the simplest version of our model has no 
uncertainty related to the fraction of the processed mate- 
rial from the ICS that pollutes the ICM. In fact, once the 
relevant measurements (ICM mass, ICM metallicity, and 
ICS mass) are made sufficiently precise, our argument 
could be turned around to place limits on the metal out- 
flow from cluster galaxies. The ICS, and the associated 
SNe, should also be a n important factor in the ene rgy 
budget of the ICM (see lBrighenti fc Mathew^ 1)19981) for 
such calcula tions relating to the SN e within the cluster 
galaxies and lDomainko et all l)2004(l for the effects of ICS 
SNe on large-scale cluster gas dynamics). Lastly, Table 
n] suggests that the factor of ~ 2 scatter in ab undances 
amon g clusters (jBaumgartner et al.l 120041 iHorner et al.l 
I2004J) is related to the balance between the mass of the 



ICM gas and that of the ICS. Wc speculate that the trend 
seen in Fc abundance versus tempe rature |Horne,r |2Q0li 
iBuote. Brighenti. fc Mathew£i20d^ may reflect an inter- 
play between the ICS fra ction, which appears t o decrease 
in less massive clusters IjFeldmeier et all [20041) . and the 
gas fraction, which varies in a complex manner as a func- 
tion of temperature (Sanderson et al. 2003). 

This paper is not a detailed exploration of the param- 
eter space currently accessible to models of ICM enrich- 
ment by the ICS, but rather a demonstration case for 
the importance of including the ICS in any accounting 
of the chemical enrichment history of clusters. As such, 
our conclusions are based on simple models that adopt 
typical values of such unknown quantities as the initial 
mass function and stellar yields. Future work should aim 
to obtain the requisite observational constraints to test 
these models in detail across a larger sample of clusters. 
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